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We investigate entanglement in the above-threshold Optical Parametric Oscillator, both theoreti- 
cally and experimentally, and discuss its potential applications to quantum information. The fluctu- 
ations measured in the subtraction of signal and idler amplitude quadratures are A 2 p_ = 0.50(1), or 
—3.01(9) dB, and in the sum of phase quadratures are A 2 g+ = 0.73(1), or —1.37(6) dB. A detailed 
experimental study of the noise behavior as a function of pump power is presented, and discrepancies 
with theory are discussed. 
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I. INTRODUCTION 

The Optical Parametric Oscillator (OPO) has been 
studied since the 196O's0 0. Already in the 1980's it 
was recognized as an important tool in quantum optics, 
for the generation of squeezed states of light 0, 0]. It 
was also recognized as a suitable system for the demon- 
stration of continuous variable (CV) entanglement, by 
Reid and Drummond, in 1988 5], where above-threshold 
operation was considered. In the early 1990's, CV en- 
tanglement was indeed demonstrated for the first time 
in an OPO, although operating below threshold @. The 
OPO has since been used in several applications in CV 
quantum information 0, @> 0- EH E3 • Entanglement in 
the above-threshold OPO, on the other hand, remained 
an experimental challenge until 2005, when it was first 
observed by Villar et al. anc ^ subsequently by two 
other groups 0, Q . 

Bipartite continuous variable entanglement can be 
demonstrated by a violation of the following inequal- 
ity, obtained independently by Duan et al. |l5| and Si- 
mon [16j: 

A 2 p_ + A 2 g+ > 2 , (1) 

where p_ = (pi - p 2 )/V2 and q + = (qi + <?2)/V2 are 
EPR-like operators constructed by combining operators 
of each subsystem. We choose pj and q~j, j € {0, 1,2}, 
as the amplitude and phase quadrature operators of the 
pump, signal and idler fields, respectively, which obey the 
commutation relations [p*j,qk] = 2iSjk- Any separable 
system must satisfy Eq. (Q: violation is an unequivocal 
signature of entanglement. 

Entanglement between the intense signal and idler 
beams generated by an above-threshold OPO can be 
physically understood as a consequence of energy con- 
servation in the parametric process. On one hand, pump 
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photons are converted into pairs of signal and idler pho- 
tons, leading to strong intensity correlations; on the other 
hand, the sum of frequencies of signal and idler photons 
is fixed to the value of pump frequency, leading to phase 
anti-correlations. The difficulty of measuring phase fluc- 
tuations was largely responsible for the long time between 
the prediction and the first observation of entanglement 
in the above-threshold OPO. The technique we used to 
measure phase fluctuations consists of reflecting each field 
off an empty optical cavity, as explained in Ref |l7j . 

The value of Eq. obtained in the first demonstra- 
tion of entanglement was 1.41(2), with squeezing ob- 
served in both EPR-like operators, A 2 p_ = 0.59(1) and 
A 2 q + = 0.82(2) Nevertheless, such a result could 
only be achieved very close to threshold, otherwise the 
phase sum A 2 q + would present excess noise, increasing 
with pump power relative to threshold a = Pq/P^. This 
strange behavior, also observed by other groups ^8|, is 
not predicted by the standard linearized OPO theory, 
for a shot noise limited pump beam. According to this 
model, entanglement should exist for all values of cr, al- 
though the degree of entanglement should decrease for 
increasing a. This presented an additional complication 
for the first demonstration of entanglement in the above- 
threshold OPO. 

In this paper, we present new improved results of en- 
tanglement in the above-threshold OPO, together with 
a theoretical and experimental study of this unexpected 
excess phase sum noise. The paper is organized as fol- 
lows. We begin by describing the linearized model for the 
OPO and its predictions for a shot noise limited pump 
beam. This model includes losses and also allows for 
nonvanishing detunings of pump, signal, and idler modes 
with respect to the OPO cavity. We then present a full- 
quantum treatment, neglecting losses and for zero de- 
tunings. Even after eliminating the linearization approx- 
imation, the theory does not predict the observed excess 
noise. The experiment is described next, and we present 
measurements of sum and difference of quadratures' fluc- 
tuations, as a function of a. The excess noise in the 
phase sum can be related to pump noise generated inside 
the OPO cavity, as we will see. We finally present our 
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currently best measurement of two-color squeezed-state 
entanglement. We conclude by mentioning applications 
of this entanglement in quantum information. 

II. THEORETICAL DESCRIPTION OF THE 
OPO 

The optical parametric oscillator consists of three 
modes of the electromagnetic field coupled by a nonlin- 
ear crystal, which is held inside an optical cavity. The 
OPO is driven by an incident pump field at frequency 
u)q. Following the usual terminology, the downconverted 
fields are called signal and idler, of frequencies u>\ and u>2, 
where, by energy conservation, loq = u)% + cl>2- We will 
treat here the case of a cavity which is triply resonant 
for wo j wi, and 0J2- Each field is damped via the cavity 
output mirror, thereby interacting with reservoir fields. 
The effective second-order nonlinearity of the crystal is 
represented by the constant x- 

Reid and Drummond investigated the correlations in 
the nondegenerate OPO (NOPO) both above [H| and 
below threshold [2(|. In the above threshold case, they 
studied the effects of phase diffusion in the signal and 
idler modes, beginning with the positive P-representation 
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where 7 and 70 are half the transmissions of the mirrors, 
7' and 7o are the total intracavity losses, fi = 7' — 7 and 
Mo = 7o — 7o ar e the spurious intracavity losses, A and 
Ao are the detunings of the OPO cavity relative to the 
fields' central frequencies, and r is the cavity roundtrip 
time. We have considered here that 71 = 72 = 7 and 
7J = 7 2 = l' ■ The parameter \ is the effective second- 
order nonlinearity. The terms 5v,j and 5vj are vacuum 
fluctuations associated to the losses from the mirrors' 
transmissions and from spurious sources, respectively. 
In the case of the intracavity pump mode, the fluctua- 
tions that come from the mirror transmission are due to 
the quantum fluctuations of the input pump laser beam, 
da n = Sp Q n + iSq n . 

Linearization consists in writing ctj (t) = e l * J (pj + 
Spj(t) + iSqj(t)) and ignoring terms that involve prod- 



equations of motion for the interacting fields [2ll |22| . 
Changing to intensity and phase variables, they were able 
to show that output quadratures could be chosen which 
exhibited fluctuations below the coherent state level and 
also Einstein-Podolsky-Rosen (EPR) type correlations. 
In the below threshold standard linearized cal- 

culation was sufficient to obtain similar correlations. In 
the limit of a rapidly decaying pump mode, Kheruntsyan 
and Petrosyan were able to calculate exactly the steady- 
state Wigner function for the NOPO, showing clearly 
the threshold behavior and the phase diffusion above this 
level of pumping j2^| . 

We begin by describing the linearized model, and then 
proceed to calculate noise spectra beyond linearization. 



A. The linearized model 

The equations describing the evolution of signal, idler, 
and pump amplitudes, aj, inside the triply resonant OPO 
cavity are given below |l7j | . They are obtained by writing 
the density operator equation of motion in the Wigner 
representation, and then searching for a set of equivalent 
Langevin equations. 
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ucts of fluctuations in the equations. Here (ctj) = p^e 1 ^ 
is each field's mean amplitude, with p\ = p 2 = p for 
equal overall intracavity losses in signal and idler, 5pj(t) 
is the amplitude fluctuation, and Sqj(t) is the phase fluc- 
tuation. Taking the average of the resulting equations 
gives information on the mean values of the fields. We 
may then separate the fluctuating part in real and imag- 
inary contributions in order to obtain the equations of 
evolution for the quadratures of the fields. Defining 
Sq± = (Sqi ± 5q2)/V% and Sp± = (Spi ± 5p2)/V% as 
the normalized sum/subtraction of signal and idler am- 
plitude and phase quadratures, we write the above equa- 
tions in terms of the EPR variables: 
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where j3 — p/po is the ratio between the intracavity am- 
plitudes of downconverted and pump fields. Noise spec- 
tra of the transmitted fields are calculated by solving the 
above equations in Fourier space. We define S p ± and S q ± 
as the noise spectra of the operators p± and q± , respec- 
tively. 

It is clear from Eq. © that the subtraction of quadra- 
tures' subspace decouples from the others, so that S q - 
and S p - depend only on the ratio of losses through the 
output cavity mirror to the total intracavity losses, and 
on the analysis frequency f2. These fluctuations do not 
depend on pump power, and are in a minimum uncer- 
tainty state, S p - x S q - — 1, if 7 = 7'. 
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FIG. 1: Prediction of the linearized theory for fluctuations in 
the sum/subtraction of field quadratures as a function of a 
for a shot noise limited pump beam. Full line: S q +; dashed 
line: S p +; line + crosses: S p -; line + circles: S q - 

On the other hand, the sum of quadratures and pump 
fields' subspaces are connected. This directly implies that 
excess noise in the pump beam degrades signal-idler en- 
tanglement, and can even destroy it ^t|- The behavior of 
the twin beams' fluctuations as functions of pump power 
relative to threshold a, for a shot noise limited pump, 
is presented in Fig. ^ The maximum squeezing of S q + 
occurs at threshold, and approaches shot noise for higher 



pump powers. 

These behaviors change in the presence of excess noise 
in the pump. In this case, both S q + and 5 P + increase 
from their values at threshold. In particular, S q + goes 
from squeezing to excess noise. The point where it crosses 
the shot noise value solely depends on the amount of 
excess phase noise present in the pump beam. For this 
reason, it was necessary to filter the pump field in the 
experiment, in order to observe entanglement. 



B. Noise spectra beyond the linearized model 

We present here a comparison between the linearized 
approach to the quantum noise in the OPO and the 
numerical integration of the quantum stochastic equa- 
tions in the positive P-representation. This will help 
us to eliminate the linearization procedure as the reason 
for the discrepancy between the theoretical prediction of 
squeezing and the experimentally observed excess phase 
noise for a > 1.2. We shall follow the procedure used 
in Ref. |24| Although exact Heisenberg equations of mo- 
tion can be found for this system, it is, at the very least, 
extremely difficult to solve nonlinear operator equations. 
We therefore develop stochastic equations of motion in 
the positive P-representation, which in principle give ac- 
cess to any normally-ordered operator expectation values 
we may wish to calculate. To find the appropriate equa- 
tions, we proceed via the master and Fokker-Planck equa- 
tions. Using the standard techniques for elimination of 
the baths |25j , we find the zero-temperature master equa- 
tion for the reduced density operator. The master ec 
tion may be mapped onto a Fokker-Planck equation 
for the positive-P pseudoprobability distribution. The 
cavity damping rates at each frequency are jj 3 = 2'yj/r, 
with 71 — 72 = 7. We further define 7 r = 70/7. In order 
to apply perturbation theory, we introduce a normalized 
coupling constant, 
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which will be a power expansion parameter. Moreover, 
it will be useful to work with the scaled quadratures 
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cc+ = gp+ , y+ = gq+, EPR variables become 

x- = gp- , y- = gq- , (5) 

in order to render the stochastic equations amenable to 
perturbation. The stochastic equations for the scaled 
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where T = 7-° i is time in units of the cavity lifetime 
for the down-converted fields. The functions £±(T) and 
£± (T) are independent Langevin forces with the following 
nonvanishing correlation functions: 



(£+(T)£ + (T')) 
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(£+(T)£+(T'))=<5(T-T'), 
(£_(T)£±(T'T)) = -<J(T - T') .(7) 



We notice the symmetry properties of the stochastic 
equations JHJ). In fact, it is easy to verify that the equa- 
tions of motion are unchanged by the transformation 
x_ <-> y + and x + <-> — y_ . Of course, all noise terms 
appearing in Eqs. © are statistically equivalent. There- 
fore, these equations should not change the symmetries 
of the initial values chosen for x+ and In order to 
provide a comparison between the linearized model and 
the full stochastic integration, we will use a perturbation 
expansion of the positive P-representation of the dynam- 
ical equations. This allows us to include quantum effects 
in a systematic fashion (S^. We first introduce a formal 
perturbation expansion in powers of the parameter g, 

oc 
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The series expansion written in this way has the prop- 
erty that the zeroth order term corresponds to the classi- 
cal field of order 1 in the unsealed quadrature, while the 
first order term is related to quantum fluctuations of or- 
der g, and the higher order terms correspond to nonlinear 



corrections to the quantum fluctuations of order g 2 and 
greater. The stochastic equations are then solved by the 
technique of matching powers of g in the corresponding 
time evolution equations. The steady state solutions Xj S 
of the zeroth order give the operation point of the OPO 
and describe its macroscopic behavior. For triply reso- 
nant operation, the expressions for the steady state are 
quite simple: 
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The first order equations are often used to predict 
squeezing in a linearized fluctuation analysis. They 
are non-classical in the sense that they can describe 
states without a positive-definite Glauber-Sudarshan P- 
distribution [23, 123 1 but correspond to a simple form of 
linear fluctuation which has a Gaussian quasi-probability 
distribution. A full quantum description of the OPO dy- 
namics can be obtained by numerical integration of the 
stochastic equations (0 , and can be compared to analyt- 
ical expressions obtained from the linearized approach. 
Taking the first order terms and using the steady state 
solutions given by Eqs. JJjJ, we can write the following 
equations for the linear quantum fluctuations, 
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The linear coupled stochastic equations obtained agree 
with Eqs. (J3J, for zero detunings and no spurious losses. 
From them, we may readily calculate the steady state av- 
erages of the first-order corrections and use that to com- 
pute the linearized fluctuations. Notice that under the 



linear approximation y_ becomes a purely diffusive vari- 
able (phase diffusion). In an experimental situation, the 
noise spectra outside the cavity are generally the quanti- 
ties of interest. We will therefore proceed to analyze the 
problem in frequency space, via Fourier decomposition 
of the fields. The first order stochastic equations may be 
rewritten in the frequency domain so that we may calcu- 
late the spectra of the squeezed and anti-squeezed field 
quadratures. The solutions for the noise of the squeezed 
operators, p_ and <?+, are: 
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where f2' = Q/j D is the analysis frequency in units of 
the cavity bandwidth. Under the limits of the linearized 
approach, the results of the noise spectra are indepen- 
dent of the phase space representation employed. There- 
fore, these results coincide with the usual ones obtained 
with the Wigner representation. The spectra given by 
Eqs. I|ll(l and i|12|) can now be compared with those 
found via stochastic integration of the full equations of 
motion JHJ in the positive P-representation. The non- 
linear spectra are calculated by Fourier transform of the 
stochastic integration, which must be performed numer- 
ically. A somewhat subtle point arises here: the non- 
linear Eqs. © have more than one possible steady-state 
solution. Thus, for a fair comparison with the linearized 
spectra, it is necessary to choose the same steady-state. 
By doing this, we verified that both predictions, in the 
above-threshold OPO, agree within a good numerical 
precision. Therefore, we conclude that possible limita- 
tions of the linearized model for dealing with the OPO 
dynamics under phase diffusion do not account for the 
experimentally observed excess noise of q+ . 



III. EXPERIMENT 

Our system is a triply resonant type-II OPO operat- 
ing above threshold. The experimental setup is depicted 
in Fig. |21 The pump beam is a diode-pumped doubled 
Nd:YAG laser (Innolight Diabolo) with 900 mW output 
power at 532 nm. A secondary output at 1064 nm is used 
for alignment purposes. Since the pump beam presents 
excess noise for frequencies as high as 20 MHz, a filter 
cavity is necessary. Our filter cavity has a bandwidth of 
2.4 MHz and assures that the pump laser is shot noise 
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FIG. 2: Sketch of the experimental setup. 



limited for analysis frequencies higher than 15 MHz (see 
Fig. |2J. We measured the laser phase noise by reflect- 
ing the beam off an empty cavity, in the same way we 
measure phase noise of the downconverted beams. The 
phase noise equals the intensity noise, except at a fre- 
quency of 12 MHz, where there is very big phase noise, 
owing to a frequency modulation inside the Diabolo laser, 
for stabilization purposes. This excess noise saturates 
our electronics and prevents measurements for analysis 
frequencies close to 12 MHz and also to its second har- 
monic, 24 MHz, as can be seen in Fig. [3J The OPO 
cavity is a linear semi-monolithic cavity composed of a 
flat input mirror, directly deposited on one face of the 
nonlinear crystal, with 93% reflectivity at 532 nm and 
high reflectivity (> 99.8%) at 1064 nm, and a spherical 
output mirror (50 mm curvature radius) with high re- 
flectivity at 532 nm (> 99.8%) and 96% reflectivity at 
1064 nm. The nonlinear crystal is a 10 mm-long Potas- 
sium Titanyl Phosphate (KTP) from Litton. Threshold 
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power is 12 mW. 
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FIG. 3: Measurement of the pump noise, as a function of the 
analysis frequency. Open circles: unfiltered laser noise; full 
circles: laser noise at the output of the filter cavity. In view 
of the large excess noise at 12 MHz and its second harmonic, 
we suppressed those frequencies from the data. 



Signal and idler beams are separated by a polarizing 
beam splitter (PBS) and sent to detection, which consists 
of a ring cavity and a photodetector (Epitaxx ETX 300) 
for each beam. Overall detection efficiency is 77 = 80(2)%. 
Both analysis cavities have bandwidths of 14 MHz, allow- 
ing for a complete conversion of phase to amplitude noise 
for analysis frequencies higher than 20 MHz. Measure- 
ments are taken at analysis frequency equal to 27 MHz. 
In order to access the same quadrature for both beams, 
the two cavities must be detuned by the same amount at 
the same time. By scanning the detunings synchronously, 
we can measure all quadratures of the twin beams. In 
particular, we can easily select the amplitude (off reso- 
nance) or phase (detuning equal to half the bandwidth) 
quadratures!^ • 

Data acquisition is carried out by a demodulating 
chain, which mixes the photocurrents from each detec- 
tor with a sinusoidal electronic reference at the analysis 
frequency and filters the resulting low frequency signal. 
The demodulated photocurrent fluctuations are sampled 
at 600 kHz repetition rate by an A/D card connected 
to a personal computer. The variances of these fluctu- 
ations are then computed taking groups of 1000 points, 
resulting in something proportional to the photocurrents' 
power spectrum at the analysis frequency. At the end, 
measured variances are normalized to the shot noise stan- 
dard quantum level (SQL). 



A. Fluctuations as a function of a 

The input pump field is guaranteed to be shot noise 
limited for frequencies above 15 MHz after being trans- 
mitted through the filter cavity. Even before being fil- 
tered, pump field is shot noise limited above 25 MHz, as 
shown in Fig. [3J Nevertheless, we observed excess noise 
in the sum of phases of signal and idler beams, prevent- 
ing the violation of the inequality given in Eq. ^ except 
for pump powers very close to threshold|l2j. 
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FIG. 4: Intensity noise of the reflected pump beam, as a func- 
tion of the detuning of the OPO cavity. The excess noise 
observed is peaked for Ao close to half the OPO cavity band- 
width. The asymmetry in the mean field signal is due to ther- 
mal bistability. The analysis frequency is 27 MHz. Circles: 
reflected pump noise; full line: reflected average intensity 

As seen in section [nj from the theoretical description 
of the OPO, excess noise in the pump beam would gener- 
ate excess noise in the phase sum of the twin beams. Yet, 
how could that be the case if we carefully measured the 
input pump to be shot noise limited? By following this 
single lead, it is natural to examine the noise properties 
of the pump beam reflected from the OPO cavity. This 
was done by scanning the OPO cavity, for crystal tem- 
peratures such that there was no parametric oscillation 
(triple resonance depends sharply on crystal temperature 
and can be easily avoided). Since the incident beam is 
shot noise limited, could there be excess noise generated 
inside the cavity containing the KTP crystal? We did in- 
deed find excess noise in the reflected pump's amplitude 
(Fig. and phase quadratures. The maximum values, 
for cr = l, were S* Q = 1.8(1) and Sf = 4.5(3). 

At present, we can still not claim to fully understand 
the origin of this excess noise. We verified, of course, 
that no such noise is generated in an empty cavity (which 
would also invalidate the measurements we perform with 
the analysis cavities for the twin beams). We also checked 
whether this effect depended on \ an d would thus be 
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directly related to the parametric process. For a polar- 
ization of the incident beam orthogonal to the usual po- 
larization, phase matching can not be fulfilled, and no 
downconversion can occur. The noise in the reflected 
beam did not show any significant dependence on the 
incident polarization, ft does, however, increase for in- 
creasing power of the incident beam. We can speculate 
that this can be a result of photon absorption by the 
crystal at 532 nm (which is at the origin of the thermal 
bistability observed in Fig.^J, with subsequent relaxation 
by spontaneous emission or non-radiative processes. This 
may give rise to an intensity-dependent refractive index, 
yielding phase and amplitude modulation at 532 nm. We 
are currently investigating these possibilities. 

As a first approximation, in order to see whether this 
would account for the behavior of A 2 p_, A 2 q_, A 2 p + , 
and A 2 (J_|_, as a function of a, we simply added ex- 
cess noise to the input pump beam in the linearized 
OPO theory. In Fig. we compare the results from 
the model, with incident S p q = 1.5 and S q o = 5.5, to 
the measured data. Signal and idler powers varied from 
0.4mW up to 5.5mW each during the experiment, corre- 
sponding to pump powers between 13mW and 26mW, or 
1.06 < er < 2.2. As expected, noises corresponding to the 
subtraction subspace, A 2 p_ and A 2 g_, are independent 
of pump power. But A 2 g + is very sensitive to <r, as well 
as A 2 p + to a smaller degree. The agreement with the 
theoretical model is surprisingly good. This is a strong 
indication that the intracavity pump excess noise is the 
main responsible for the excess noise in A 2 q + . 



B. Two-color entanglement 

The sum of phases' noise is squeezed very close to 
threshold, and squeezing is degraded with increasing 
pump power. A 2 q + crosses the shot noise level approx- 
imately at a = 1.20, from squeezing to anti-squeezing, 
although only below a = 1.15 can squeezing be observed 
with certainty. 

Fig.|S|shows the recorded noise in sum and subtraction 
of photocurrent fluctuations of signal and idler beams as 
functions of analysis cavities' detuning, for a = 1.06. Off 
resonance, quantum correlations are observed in the sub- 
traction of amplitudes, A 2 p_ = 0.50(1), or -3.01(9) dB. 
For analysis cavities' detuning equal to half the band- 
width, squeezing is present in the sum of phases, A 2 q + = 
0.73(1), or -1.37(6) dB. The Duan et al. and Simon cri- 
terion, Eq.QJ, is then clearly violated, 
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A 2 q+ = 1.23(2) < 2 



(13) 



attesting the entanglement. This value, together with 
the one reported by Jing et al. 01 > is the lowest achieved 
for twin beams produced by an above-threshold OPO. 

We also point out that, in this experiment, the twin 
beams have very different frequencies (wavelengths dif- 
fer by ss 1 nm), an unusual situation. Such two-color 
entanglement can be very interesting for the transfer of 
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FIG. 5: Noise behavior as a function of a. In part (a), we 
present the predictions of the linearized model, for an input 
pump beam with S p o = 1.5 and S q o = 5.5; dashed line: 5^+; 
full line + open circles: S q -; full line: S q +; full line + crosses: 
S p -\ SQL = 1.0 is indicated by a dashed line. In part (b), 
experimental results are shown for a ranging from 1.06 to 2.2. 
full circles: S p+ ; triangles: S q -; open circles: S q+ ; squares: 
S p _; SQL = 1.0 is indicated by a dashed line. 



quantum information between different parts of the elec- 
tromagnetic spectrum. 



IV. CONCLUSION 

We presented a theoretical and experimental investi- 
gation of phase noise and entanglement in the above- 
threshold OPO. Excess noise in the phase sum of the 
twin beams was measured as a function of pump power 
relative to threshold and we found that it decreases as 
pump power is lowered. We finally discovered that excess 
pump noise is generated inside the OPO cavity contain- 
ing the nonlinear crystal, even for a shot noise limited 
pump beam and without parametric oscillation. The ul- 
timate physical origin of this phenomenon still requires 
further investigation. Another important question to ad- 
dress is how one can eliminate this effect. Su et al.fl^ 
were able to observe entanglement for a of the order of 
two. The difference between their setup and others is a 
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FIG. 6: Sum (full circles) and difference (open circles) of 
quadratures' noise, measured as a function of the analysis 
cavities' detuning. Squeezed-state entanglement can be di- 
rectly observed, with A 2 p_ = 0.50(1), or -3.01(9) dB, and 
A 2 q+ = 0.73(1), or -1.37(6) dB 

lower cavity finesse for the pump field. If the assumption 
of an intensity dependent index of refraction is correct, 
this makes sense. For a lower finesse, phase shifts accu- 
mulated inside the cavity should be smaller, hence the 
excess noise generated should also be smaller. 

In spite of these unexpected phenomena, two-color en- 



tanglement was measured in the above-threshold OPO. 
There are interesting avenues to pursue for applications 
in quantum information. First of all, we should mention 
that the strongest squeezing measured to date, —9.7 dB, 
was generated in an above-threshold OPOhJ. Thus, 
entanglement in the above-threshold OPO may be the 
strongest ever achieved for continuous variables. The 
bright twin beams can have very different frequencies, 
and one can envisage CV quantum teleportation^ to 
transfer quantum information from one frequency to an- 
other (in other words, to "tune" quantum information). 
For example, this system could be used to communi- 
cate quantum information between quantum memories 
or quantum computers based on "hardwares" which have 
different resonance frequencies. Finally, a quantum key 
distribution protocol proposed by Silberhorn et aL[3l| 
can be readily implemented, with the advantage that the 
measurement with analysis cavities does not require send- 
ing a local oscillator together with the quantum channel 
to the distant receiver. 

The above-threshold OPO, which was the first system 
proposed to observe continuous variable entanglement, 
has finally been added to the optical quantum informa- 
tion toolbox. We expect new and exciting applications 
to come in the near future. 
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